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Abstract 

This paper describes a new 2D model for the photospheric evolution of the 
magnetic carpet. It is the first in a series of papers working towards construct- 
ing a realistic 3D non-potential model for the interaction of small-scale solar 
magnetic fields. In the model, the basic evolution of the magnetic elements 
is governed by a supergranular flow profile. In addition, magnetic elements 
may evolve through the processes of emergence, cancellation, coalescence and 
fragmentation. Model parameters for the emergence of bipoles are based upon 
the results of observational studies. Using this model, several simulations are 
considered, where the range of flux with which bipoles may emerge is varied. In 
all cases the model quickly reaches a steady state where the rates of emergence 
and cancellation balance. Analysis of the resulting magnetic field shows that we 
reproduce observed quantities such as the flux distribution, mean field, cancella- 
tion rates, photospheric recycle time and a magnetic network. As expected, the 
simulation matches observations more closely when a larger, and consequently 
more realistic, range of emerging flux values is allowed (4 x 10^^ — 10^^ Mx). The 
model best reproduces the current observed properties of the magnetic carpet 
when wc take the minimum absolute flux for emerging bipoles to be 4 x 10^^ 
Mx. In future, this 2D model will be used as an evolving photospheric boundary 
condition for 3D non-potential modelling. 

Keywords: Sun: magnetic fields - Sun: magnetic carpet 



1. Introduction 

The solar magnetic carpet is a salt and pepper mix of small-scale opposite 
polarity magnetic features. These features are continually evolving in response 
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to the photospheric flow pattern of convective cells known as supergranules. 
These cells range in diameter from 10 000 km to 50 000 km, with an average 



of between 13 000 and 18 000 km (Hagenaar et al., 19971. Their flow pattern is 



described by an upflow at cell centre, followed by radial outflow at roughly 0.5 km 
s~^ ( [Simon and Leighton, 1964} [Simon et al, 2001] |Paniveni et al., 20041, and 



downflow on the order of 0.1 km s~ at the cell edges ( Wang and Zirin, 1989 1. 

Due to the constant motion of the underlying granules and supergranules, 
magnetic features frequently interact with one another. New magnetic flux is 
continually emerging from below the surface in the form of opposite polarity 
pairs ( [Harvey and Martin, 1973 1 . In addition, other features are formed through 



the coalescence of same-polarity features; large features may break apart into 
smaller features; and features may shrink or disappear through the cancellation 
of opposite polarity flux. Over the last 10 years our understanding of the solar 
magnetic carpet has greatly increased due to the space missions of SOHO/MDI 
and Hinode /SOT. Through these missions, many observational studies of the 



properties of the magnetic carpet have been carried out (Schrijver et al., 1997 
Hagenaar, 2001|[Parnell, 2002| [Hagenaar et al., 2003|[DeForest et al, 2007|[McIntosh et al., 2007|| 



de Wijn et al, 2 008"Ha.gcna ar et al, 2008[[Lamb et al., 2008[[Lites, 2009[[Parnell et al, 2009[ | 



Lamb et al., 2010; Thornton and Parneh, 20111 



Three main types of small-scale magnetic features have been classified; these 
are ephemeral regions, network features and internetwork features. Ephemeral 
regions are newly emerging bipolar pairs that appear within the supergranular 
cells. They have an average lifetime of 4.4 h and typical diameters of 3 000- 
5 000 km ( [Harvey and Martin, 1973 ). For the first 30 mins after appearing, the 



two opposite polarities are found to separate from one another at a velocity 
of around 4.5-5 km s~^ ( [Harvey, 1993[ [Title, 2000[ ). Later, they slow to a ve- 
locity that is on the order of the underlying supergranular flow, « 0.4 km 
s~^ ( [Hagenaar et al, 2003p . [Schrijver et al.\ (|1997p give the average flux of an 
ephemeral region to be 1.3x10^^ Mx (Maxwell), while IChae et a/.l([MJT)) find a 
slightly higher value of 2.8x10^9 Mx. 

Network features are larger features, with diameters in the range 1 000-10000 
km and fluxes of 10^^-10^^ Mx, that are typically found at sites of strong down- 
flow, i.e. the edges of the supergranular cells ( [Martin, 1988[ ). IZirinI (jl985p find 
them to be slow moving, with an average velocity of just 0.06 km s~^. They do 
not emerge as network concentrations; rather they are produced from the resid- 
uals of other magnetic flux features. Around 90% of their flux originates from 
ephemeral regions, with the remaining 10% arising from internetwork features 
( [Martin, 1990| ). 

The flnal type of small-scale flux features are internetwork features; they 
also emerge as bipolar pairs within supergranular cells. Their mean diameter is 
around 2 000 km and their fluxes lie in the range 10^^-2 x 10^^ Mx ( [Wang et al, 1995[ ).| 



De Wijn et al. (|2008p find a rms velocity of 1.57±0.08 km s~^, and an average 
lifetime of just 10 mins for internetwork features. IZhou et aD (|2010p deduce an 
even smaller average lifespan of just 2.9±2.0 mins. 

Previous theoretical models for the evolution of the magnetic carpet photo- 
sphere include [Schrijver et al.\ (|1997p , [van Ballegooijen et al.\ (|1998p , ISimon et al. ^ 
(pngT|) . lPaFSen(PDnT|) and [Cranmer and van Ballegooijen] | (|2{mi)) . [ScE-ijver et al.\ 



SOLA: km_paperl.tex; 5 August 2011; 0:09; p. 2 



Simulation of Synthetic Magnctograms 



(|1997p describe flux emergence, cancellation, coalescence and fragmentation througli§ 
equations of magneto-chemistry. Other authors represent magnetic flux fea- 
tures as collections of corks that are passively advected by photospheric flows 



van Ballegooijen et al., 1998 Simon et ai, 20011, or as a continuous distribu- 



tion ( Parnell, 2001 Cranmer and van Ballegooijen, 2010 ). 

This paper describes a new model that aims to reproduce the main observa- 
tional properties of this small-scale photospheric field. Many of the parameters 
built into it are taken from studies of observational data (such as Hinode/ SOT 
or SOHO/MDI magnetograms) . The model also includes a supcrgranular flow 
pattern ( [Simon and Leighton, 1964P , as well as random small-scale motions to 
represent the effect of granulation. 

The model reproduces observed quantities, such as the power law flux distribu- 
tion obtained bv lParnell et aLI(|2009p . Other objectives include reproducing the 
rapid photospheric recycle time on the order of 1-2 h deduced by |Hagenaar et aT] 
(|2008p : equality of the rates of emergence and cancellation of flux; and visually, 
the appearance of a magnetic network along the boundaries of the supergranule 
cells. 

Our long term goal is to use the simulated magnetograms produced by this 
model as lower boundary data in our 3D non-linear force- free coronal field model 
( [Mackay et ai, 2011[ ). This model applies a magnetofrictional method to pro- 
duce a non-potential coronal magnetic field in response to photospheric motions 
dYeates et ai, 2008[ ). The method wfll be described in Paper 2 ( [Meyer et ai, 2011bp ,| 
in which we simulate the non-linear force-free evolution of the coronal field for 
a scries of simple small-scale interactions. 

The paper is outlined as follows: Section [5] describes the technical aspects of 
the construction of the 2D magnetic carpet model and how each of the four flux 
evolution processes are specifled. The results of the simulations are presented 
in Section [3] where we show that the model reproduces many observational 
parameters. Section |3| gives the discussion, conclusions and an outline of future 
work we are undertaking. 



2. Model 

The construction of the magnetic carpet model is now described. Whilst ob- 
servational studies use the term 'feature' to describe small-scale magnetic flux 
concentrations, we will refer instead to magnetic 'elements'. This is to distin- 
guish between a 'magnetic element' in the model, which we choose to have a 
specific mathematical form. In contrast, an observational magnetic feature as 
would appear in a magnetogram may be composed of several of our magnetic 
elements. Each magnetic element within the model is treated individually as 
a unique discrete element. The sum of all elements then produces a synthetic 
magnetogram. This approach differs from other methods in several ways. In cork 
models, each cork represents a single intense flux tube, whereas in our model 
we represent our magnetic elements with Gaussian peaks, allowing features to 
form that are composed of many peaks and troughs in magnetic field strength. 
To avoid undesirable numerical effects such as numerical diffusion and pilc-up 
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Figure 1. Diameter of a magnetic element versus absolute flux. 

at cancellation points, we move the centres of magnetic elements rather than 
advect their Gaussian profiles. Such a treatment also allows us to easily keep 
track of the number of elements and exactly which elements are involved in each 
of the four processes of emergence, cancellation, coalescence and fragmentation 
at any time. We first describe the mathematical form of the magnetic elements 
which produce the synthetic magnetograms. Following this, examples of the 
magnetograms produced over a 250 h period are shown. Finally, we discuss the 
rules that govern how the magnetic elements evolve. 

2.1. Synthetic Magnetograms 

For each discrete magnetic element we assume that the z-component of the 
element's magnetic field has a Gaussian profile. 



where Bq is the peak magnetic field strength, tq is the Gaussian half-width and 
r is the distance from the centre of the Gaussian. The total flux, $, of each 
element is: 



J A 

We specify the Gaussian half-width of each magnetic element be Vq = (i((/))/4, 
where d{(j)) is the diameter of the magnetic element and (j) is its absolute flux, 
(/) = I $ I . The diameter is given by 



(1) 




= TO0logio('/') 



(3) 



where the parameters and arc given by 



dmax ^ 'j^min 



and C0 = djjjiji 



m<t, logio(0small)- 



7710 = 



logio('/'large) -loglo('/' small ) 
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Figure 2. Synthetic magnctograms for a simulation with a flux emergence range of 8 X 10^'' 
- 1 X 10^® Mx. Both are taken at t = 20 h. (a) Red contours represent positive magnetic field, 
blue represents negative magnetic field, where 10 contour levels are shown for each polarity 
with an absolute peak value of 70 G (Gauss), (b) The same as image (a), with contours filled 
in to simulate a magnetogram. 



We let d^jy[Yi = 1 Mm, dmax = 6 Mm, ^gmall = ^^^^ ^^^'^ ^large = ^^^^ 
Mx so that the majority of magnetic elements within the simulation are con- 
fined to the range d G [1,6] Mm. These values agree with observed diameters 
and fluxes for magnetic carpet features such as small network features and in- 
ternetwork features ( [Harvey and Martin, 19731 [Martin, 1988|[Wang et al, 1995[ 
[Zhou et al., 2010[ ) . A plot of diameter versus flux is shown in Figure [Tj 

Equation <^ may be rearranged to give an expression for the peak magnetic 
field strength of a single magnetic element. 



16$ 



The contribution of the jth magnetic element to the normal component of the 
magnetic field is then: 

16r2 
d(^' 



Bz,j = Bo,j exp <{ \ , (5) 



r = {x ~ Xj)- + {y - y^) , 

where {x,y) is an arbitrary position, and {xj,yj) is the position of the centre 
of the magnetic flux element. We sum up the contribution from every magnetic 
element to give B^, 

N 
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Figure [5] shows contour plots of Bz for a simulation in which newly emerging 
bipoles have a total flux in the range 8 x 10^^ -fx 10^^ Mx. This simulation 
is described in more detail in Section [31 A black and white image of the same 
region is shown to simulate a magnetogram. 

Throughout the simulation, the motion of the magnetic elements is deter- 
mined by the underlying supergranular flow, which is described next. In addition, 
details of the flux emergence and interaction processes are given in Sections 12. 3L 
12.41 and 12.51 Readers not wishing to read the technical details of how these 
processes arc implemented may jump to Section [3] 



2.2. Steady Flow Profile 



We construct a supergranular velocity profile for our model that is similar to 
that of ISimon et al.\ (j2001|) . except that the diameter of each cell varies from 
supergranule to supergranulc. For simplicity our flow profile is currently steady 
throughout the simulation, although in reality, supergranules evolve in time. 
Various authors have estimated lifetimes of supergranules to be anywhere from 10 
hours to 2 days, depending on the technique used ( [Rieutord and Rincon, 20"T0l ). 
However, wc will show that our steady profile does not lead to an unphysical 
buildup of magnetic flux at the cell boundaries. The flow profile of a single 
supergranule is given by 

■yij = ^o-Rcxp |-^|, (7) 

where vr is the velocity from cell centre in the x-y plane. R is the distance 
from cell centre, and and Rq are normalised values that determine the flow 
strength and radius of the supergranule. We set Rq = 0.173 for all supergranules, 
which corresponds to 8.7 Mm in our simulations. Aq is a number between and 1 
that determines the strength of each cell. We choose the positions (ccc, yc) of the 
centres of n supergranules in the simulation. To introduce the influence of cells 
outside the computed domain, we translate the positions of these n cells eight 
times to surround the original pattern, as illustrated by Figure [Sj^a). The original 
domain is shown in black, and spans the range [xjjjjjj, Xmax] x [?/min' J/max] , where 
we choose x^^^ = 2/jjiin = Mm and Xmax = 2/max = 50 Mm. Translation of the 
cells means that the supergranular flow matches through the side boundaries, 
which are periodic. 

To produce the steady fiow pattern, the velocity at any point is found by 
summing over the contributions of 9n supergranular cells. This includes the n 
original supergranules which describe the central region of Figure [3l^a), outlined 
in black, and the eight sets of translations. Once the contribution of all cells 
has been calculated, the resulting velocity profile is scaled so that the maxi- 
mum value of vr is 0.5 km s~^ (Simon and Leighton, 1964 Paniveni et ai, 2004 



[Rieutord and Rincon, 2010] ). 

The flow profile chosen for the simulation illustrated in Figure [3| is shown in 
Figure [3]Jb). This flow profile is used for all of the simulations described in this 
paper. The positions of the cell centres were selected to produce an irregular 
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Figure 3. (a) The positions of the supcrgranule cell centres are translated eight times to 
surround the originals. This simulates the influence of outside flows on the domain, (b) A sim- 
ulated supergranular flow profile in which seven cells have been specified in the computational 



pattern, as would be seen in real observations of the solar surface. One can also 
see flows contributing from the supergranules through the boundaries, as a result 
of the translation previously described. 

Contributions from granular motions are also added to each magnetic ele- 
ment's velocity at each time step. The granular velocity contribution is randomly 
chosen to be between and 0.1 km s~^ if an clement is within 0.75i?sg of the 
centre of the supergranular cell, and between and 0.2 km s^^ if it is further 
out than this, where i?sg is the radius of the supergranule. The direction of 
this velocity is also chosen as random. A slightly higher value for the granular 
velocity contribution is chosen if the magnetic element is near the boundary of 
the supergranule, since the contribution from the supergranular flow profile is 
small at these locations. This prevents elements from becoming stationary once 
they reach the network between supergranules. 

One time step within the model is equal to 1 min. This is an arbitrary choice 
within the model and may be chosen to be higher or lower. We have chosen 
1 min as the current cadence of magnetogram data from instruments such as 
SOHO/MDI and SDO/HMI are of a similar size. Fuh disc MDI magnctograms 
arc typically of cadence 1 min at best ( [Scherrer et al., 1995[), while HMI line-of- 
sight magnetograms are of cadence 45 s (http:/ /hmi. stanford.edu/hmi-specs. html 1. 1 
Many authors studying the observational properties of the solar magnetic carpet 
have used magnetogram data sets of cadence roughly 1 min, occasionally aver- 
aged over a longer time period to reduce noise, for example |Hagenaar| (|200ip , 
IParnelll ([M?^ . IZhou et al.\ (PHTU)) and [Thornton and Parnelll (POTTj) . As higher 
cadences become available for observational data, we may choose a smaller time 
step in our model. The simulations in Section |3] are each run for 250 h, which 
corresponds to 15 000 time steps. 



SOLA; km_paperl.tex; 5 August 2011; 0:09; p. 7 



Meyer et al. 




W"^ 10" 10 

Flux (Mx) 



Figure 4. Log-log plot of the frequency of emergence against flux emerged in a 5 h long set of 
Hinode/SOT magnetograms. Results were obtained using a tracked bipolar (TB) and tracked 
cluster (TC) method applied to clumping (subscript 'c') and downhill (subscript 'd') data. See 
Thornton and Parnell (2010) for definitions of tracking methods used. The dashed line is a 
power-law fit as described by Equation with no = 1.77 X 10^^* cm^^ day^-"^ and a = 2.74. 



In the next three sections, we discuss the methods used to implement the 
processes of emergence, fragmentation, cancellation and coalescence. 

2.3. Emergence 

Each emerging bipole is made up of two separate magnetic elements that are of 
equal flux and opposite polarity. We use the term 'bipole' only when referring 
to a newly emerging pair of magnetic elements, at all other times we consider 
single magnetic elements. 

To simplify computations and prevent infinitesimally small flux elements from 
arising, we set the minimum unit of flux that a single polarity may have to be 
equal to (po = 10^^ Mx, and define all magnetic elements to have flux that is an 
integer multiple of 0o ( jParnell, 2001[ ). 

A bipole of absolute flux 0bp consists of two magnetic elements of equal 
absolute flux, = ^bp/S, but opposite polarity. The total signed flux of the 
bipole is therefore 0. 

2.3.1. Parameters for Newly Emerging Bipoles 

Flux emergence within the simulation is determined by the probability distribu- 
tion for emerging bipoles from [Thornton and Parnelll (|201ip . This is determined 
from Hinode/SOT high resolution magnetograms and is in the form of a power 
law: 
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where uq = 1.77 x 10^^'' cm^^ day^^ and a = 2.74. The values for uq and a 
come from a feature trackmg study of Hinode /SOT data that is described by 
IThornton and Parnelll (|201ip . Figure H] shows a log- log plot of their results for 
the frequency of emergence versus flux emerged. TB and TC indicate emergence 
detected by a tracked bipolar and a tracked cluster method respectively. These 
methods are explained fully in the original paper. The subscripts 'd' and 'c' 
stand for the downhill and clumping methods of magnetic feature identification, 
these are described by IDeForest et al.\ (|2007p . 

The quantity N{(j)y^p)d(j)\^p is the total number of bipoles that emerge with 
total absolute flux in the range [4>hp, '/'bp + c^^bp], where c?0bp is very small. The 
flux emergence rate for bipoles with flux in the range [(j)a, ^b] in Mx cm~^ day~^ 
may be computed as: 



jV(0bp)<Pbp«(pbp = 7: 

2 - a 



"bp 

00 



(9) 



Correspondingly, the number of bipoles emerging in the range 0b] in units 
of cm~^ day~^ is then 



iV(0bp)rf0bp 



00 V 00 / I - a 



"bp 

00 



(10) 



Let 



and 0max be the minimum and maximum flux allowed for emerging 



bipoles in our simulation. They must both be integer multiples of 20o, as each 
individual polarity will then have an absolute flux that is an integer multiple 
of 00. In addition to this, we must have (t>^[j^ > 20o- For our simulations, we 
choose to emerge only discrete values of flux from 0ijiin to 0max in steps of size 
d0 = 200. If m = (0max — 0min)/20o then the set of emerging flux values of 
bipoles is 

{0fc = 0min + 2fc0o, 0, 1, m}. 

For each discrete value 0^, we integrate over the range [0^ — 0o,0fc + 0o] to 
approximate the number of bipoles with that absolute flux that will emerge 
during the simulation: 



4>k+4>a 



00 



"bp 

00 



"bp 



X Ax 



(11) 



where A is the area of the domain in cm^ and D is the number of days of the 
simulation. We cannot allow a non-integer number of bipoles to emerge, so iV^ 
must be converted to an integer, Nk, by rounding up or down randomly. The 
N = X^fcLo bipoles are then randomly assigned a value from to imax = 
15 000 using a uniform distribution, which will be the time step in which they 
emerge. This means that a random number of bipoles emerge each time step, 
and also a random total quantity of flux. 

In addition to flux, diameter and sign, several other parameters must be 
chosen for each newly emerging bipole. For each bipole, a random uniform integer 



SOLA; km_paperl.tex; 5 August 2011; 0:09; p. 9 



Meyer et al. 




(d) (e) (f) 

Figure 5. A sequence of still images (left to right) showing a newly emerging bipole (elements 
2 and 3). Red contours represent positive magnetic field, blue contours represent negative 
magnetic field. 20 contour levels are shown for each polarity, with an absolute peak value of 70 
G. The region is 10 Mm X 10 Mm in area. The images are taken 4 mins apart, and elements 
of interest have been numbered from 1 to 3, where elements 2 and 3 illustrate emergence. This 
event may be seen in the accompanying movie, 'mag2_em.mpg'. 

between and n — 1 is chosen. This is the index of the supergranule cell that 
it will emerge in. Its location within the cell is then randomly chosen. Since 
ephemeral regions have been observed to emerge with a preference towards the 
edge of a supergranule QWang, 1988[ ) , we also build this into our model. Wc allow 
emergence to occur in the range [0.5i?sg, 0.75i?sg], where i?sg is the radius of the 
supergranule. The final parameter that must be chosen for an emerging bipole 
is its tilt angle, 9. This is the angle of the axis along which the two polarities 
of the bipole will separate from one another and it is simply a random uniform 
number between and 27r radians. The separation velocities of newly emerging 
magnetic elements are discussed below. 

2.3.2. Appearance of Newly Emerging Bipoles 

At each time step, newly emerging bipoles are added into the simulation. Each 
magnetic element undergoes the specified emergence process outlined below until 
it has travelled e^^^ x d{<f)) ~ l.5d{4>) from its initial position. Cj-j^d 1-5 is a 
constant that defines the separation that a magnetic element must reach from its 
emergence point before supergranular fiows or other processes may take over its 
evolution. The element's velocity v^ep depends upon how long it has been emerg- 
ing for and its direction is given by its tilt angle 9. Observational studies show 
that the two polarities of a newly emerging bipole initially separate at several 
km s^^, later slowing to velocities that are of a similar order to those of the un- 
derlying supergranular flows ( |Harvey, 1993| [Title, 2000| [Hagenaar et al., 2003[ ). 
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Within our simulations, two elements will initially separate with a velocity of 3 
km s^^, slowing to 1.0 — 1.3 km s^^ after 30 mins and later to 0.5 km s^^ which 
is on the order of the underlying supergranular flow. Initially, the positive and 
negative magnetic elements within the bipole move in opposite directions along 
the axis of their tilt angle. 

FigureO (and the accompanying movie, 'mag2_em.mpg') shows a sequence of 
images taken from a small region of the synthetic magnctograms produced by one 
of the simulations. The region is 10 x 10 Mm^ in area, and the images are taken 4 
mins apart. One can see that between images (a) and (b), a new bipole begins to 
emerge in the centre of the box, the positive and negative polarities are marked 
'2' and '3' respectively. The two polarities grow in flux as they move apart. 
Towards the end of the sequence, magnetic element 3 begins to interact with a 
pre-existing magnetic clement, indicated by '1'. Notice also that in addition to 
the two polarities moving in opposite directions to one another, there is a slight 
drift of the bipole towards the upper left due to the underlying supergranular 
flow. 

2.4. Fragmentation 

Fragmentation within our model is based upon the process described bv lParnelll 
(|200ip . which depends upon both the flux of the magnetic element and time. 
Every element is checked for fragmentation at each time step. The fragmentation 
rate is an input parameter for the simulation. We assume that every element 
of sufficient flux will fragment within Tf = l/i?f s. IParnelH (|200ip suggests that 
a fragmentation rate of i?f > 1.2 x 10^"* s^^ is required to reproduce the correct 
absolute flux density and flux distribution, so we take i?f = 1.5 x lO"* s'^. This 
means that a magnetic element of sufhcicntly large flux will fragment roughly 
once every 1/ Rf = 6667 s w 1 h 50 mins. For an element j to have sufficient flux 
to fragment, it must satisfy the inequality 

</.,>^(l-^), (12) 

where g is a random number such that kf < q < 1, and < fcf < 1. fcf is fixed at 
the start of the simulation, but a random q is chosen for every magnetic element 
every time we check whether it will fragment. IParnelll (I200ip take a value of ip 
that is at the large end of their range of expected flux values, ^jj = 8 x 10^^ Mx, 
and /cf = 0.75. Without observational evidence to suggest differently, we take 
the same values for our simulations. With these parameters, the largest value 
that the right hand side of Equation can take is 

Vfl ^) = 8x1018(^1-^^^ = 2 X 10^8 Mx. 

V max(g)y V 1-0/ 

Similarly, the minimum value that the right hand side of Equation ()12p can 
take is 0. However, elements are not allowed to fragment if they are of unit flux, 
(j)o- This means that all magnetic elements of absolute flux > 2(j)Q may fragment, 
but elements of greater flux have a higher probability of fragmentation. 
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Time dependence is built into the fragmentation process as follows. If 0j 
satisfies Equation p^ . we then choose a random s such that < s < 1. Element 
j will only fragment during the current time step if 



Tj is the 'age' of magnetic element j in s. An element's age is reset to every time 
it fragments, coalesces or cancels with another element. Clearly once Tj > 
the right hand side will be greater than 1, and the inequality will be satisfied. 
Therefore, all magnetic elements of flux greater than 2 x 10^® Mx are guaranteed 
to fragment within T = Tf/60 time steps, unless some other process takes over 
their evolution before then. 

2.4-1- New Elements Arising From Fragmentation 

Within our simulations, elements are allowed to split into just two new elements 
at a time. In reality, most fragmenting magnetic features only split into two 
( [Thornton, 2011[ ). Once it has been determined that element j will fragment, 
a new element k is introduced to represent the element resulting from the 
fragmentation. The flux and diameter of j and k must now be recomputed. 
As in lParnein (|200ip . the original flux (pj is split into two new fluxes, cf>i and 4)2- 
If p is a random number between 0.55 and 0.95, then 

01 = P(f>] ■ 

Since all fluxes within the simulation must be integer multiples of we round 
<j>i/(f>o to the nearest integer. We then take 

02 = 4>j - 4>i- 

Elements of flux (f>Q arc not allowed to fragment. If (f>j — 2(j)Q, then 01 = 02 = 0o, 
and if 0^ = 30o, then (pi = 20o and 02 = 0o- We may now set 0^ = 0i and 
0fe = 02; then compute their new diameters d{(f>j) and d{(j)k)- 

2-4-2- Motion of Fragmenting Elements 

An element's behaviour after it has fragmented depends upon the process that 
the original element was undergoing before the fragmentation occurred. If it was 
undergoing emergence, fragmenting from yet another element, or simply being 
advected by supergranular flows, then its velocity returns to being determined 
by this proces^. In the case of emergence, the two new elements produced 
by fragmentation are also treated as emerging. Their velocities therefore still 
decrease with time, though the directions of motion of the two new elements j 
and k are different. If 9 is the direction of motion of the original element, then 

01 = 9 + 0.37r«7 and 92^9- 0.37rg, 



^If the underlying supergranular velocity is very small (< 0.1 km s the new elements may 
be given a 'push' of 0.1 — 0.2 km s^^ to help them separate. 
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Figure 6. A sequence of still images (left to right) that features fragmentation of magnetic 
elements. Red contours represent positive magnetic field, blue contours represent negative 
magnetic field. 15 contour levels are shown for each polarity, with an absolute peak value of 70 
G. The region is 10 Mm X 10 Mm in area. The images are taken 2 mins apart, and elements of 
interest have been numbered from 1 to 4. These interactions may be seen in the accompanying 
movie, 'mag3_frag.mpg'. 



where < g < 1 is randomly chosen for each individual fragmentation ( Parnell, 2001[ ).[ 
Therefore the two magnetic elements move in the same direction with slightly 
offset motions. 

If the element determined to fragment was undergoing cancellation or coa- 
lescence with another element I before it fragmented, the treatment is slightly 
different. The fluxes and diameters of the two new elements j and k are com- 
puted. The new element j then continues to cancel or coalesce with I, while k 
moves off in a different direction at the velocity of the underlying supergranular 
The new direction of motion for k is given by 

6*2 = 6* ± 0.67rq, 

with < q < I random. 

The parameter f^-^^ defines the separation that a magnetic element will reach 
from its point of fragmenting. An element j will continue to move away from its 
fragmentation site until it has reached a distance of /rad'^(<^j) = 1.5c?(0j). 

Figure [6] (and the accompanying movie, 'mag3_frag.mpg') shows examples 
of fragmentation occurring within one of the simulations. Element 1 splits to 
form a new element, 3, element 2 then begins to fragment to form element 4. 
Fragmentation may also be seen in the sequence in Figure [5] of the following 
section, where element 2 fragments into a large and a small element. 

2.5. Coalescence and Cancellation 

Every time step, each element that is not currently undergoing the emergence 
process is checked to see if it will cancel or coalesce with another element. In 
the current model, an element j may only cancel or coalesce with one single 
element at any given time. First we determine how many magnetic elements are 
within interaction range of j (including those through the periodic boundaries). 



^Again, the element may be given a 'push' if the underlying supergranular velocity is very 
small. 
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(a) (b) (c) (d) (e) 

Figure 7. A cancellation between two magnetic elements of equal flux and opposite polarity. 
The bipole shrinks as the elements move together and their Gaussian profiles overlap. By the 
last frame the magnetic elements overlap completely and no contours of flux remain. 

An element k is defined to be within interaction range if d^ep, the separation 
distance, satisfies 

4ep(j, k) < Cj.a(j(d((^j) + d{(j)k)), 

where 

dsep{j,k) = \{xj,yj) - {xk,yk)\ and Cj-ad = 0-5- 

The constant Cj-^d defines the interaction range between two magnetic elements. 
It is important that Cj-^j < 0.75 to prevent elements that have just sepa- 
rated from one another due to emergence from immediately cancelling with 
one another again. There are also several other conditions for the cancelling or 
coalescing of j with some element fc: 

• k must not currently be undergoing emergence, 

• k must not currently be cancelling or coalescing with another element, 

• If j and k are in the process of fragmenting from one another, then they 
are not allowed to immediately coalesce. If they are still within range after 
the fragmentation process has ended, they may then coalesce together (see 
Section HH). 

If these conditions are all satisfied by more than one element within range of j, 
then j will cancel or coalesce with the closest element. 

Once it has been determined that two elements will cancel or coalesce, they 
move together at a constant velocity of Vq = 1-0 km s~^ until their centres meet. 
Their direction of motion is along the axis defined by their centres. 

Since each magnetic element is given a Gaussian profile, their profiles overlap 
as they move towards one another. In the case of cancellation, this results in 
both magnetic elements shrinking as they move towards one another. If they are 
of equal flux and opposite polarity, they will completely cancel one another out. 
This is illustrated in Figure[71 Within the magnetic carpet model, both magnetic 
elements are removed at this stage. In the case of partial cancellation or coales- 
cence, one of the elements is removed when the peaks of their Gaussian profiles 
meet. The remaining element's flux is then updated to be the difference between 
the two original flux values. This new element's motion is now determined by 
supergranular flows until another process takes over. 
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(e) 



(f) 



(g) 



(h) 



Figure 8. A sequence of still images (left to right) showing magnetic elements coalescing and 
cancelling. Red contours represent positive magnetic field, blue contours represent negative 
magnetic field. 15 contour levels are shown for each polarity, with an absolute peak value of 70 
G. The region is 10 Mm X 10 Mm in area. The images are taken 4 mins apart, and elements of 
interest have been numbered from 1 to 6. These interactions may be seen in the accompanying 
movie, 'magl_canc.mpg'. 



The processes of cancellation and coalescence occur frequently within the 
model, and can easily be detected by eye in the synthetic magnetograms. Fig- 
ure[5](and the accompanying movie, 'magl_canc.mpg') shows a sequence of eight 
images taken from one of the simulations, in which examples of both cancellation 
and coalescence of magnetic elements can be seen. We follow the evolution of a 
negative polarity element 1, and two positive elements, 2 and 3. In frame (d) an 
emergence occurs in the bottom left (5 and 6) and element 2 begins to fragment, 
creating a new element, 4. As the sequence progresses, a cancellation occurs 
between 1 and 2, and elements 3 and 4 coalesce. By the end of the sequence, 2 
and 4 have disappeared, and 3 has begun cancelling with 5. 

2.6. Summary of Photospheric Evolution 

After the parameters for newly emerging bipoles have been determined, we com- 
pute the magnetic field on the photosphere for each time step. At each time step, 
we first loop over all magnetic elements and determine their current velocities. 
In addition to the four flux evolution processes described in Sections 12.31 12.41 
and 12.51 we add a contribution from the underlying supergranular flow to each 
element's velocity. This prevents the motion of elements from being too linear, 
particularly when undergoing emergence. It seems likely that supergranular flows 
would continually have an effect on magnetic elements on the real solar surface. 
We also add a random motion to every magnetic element that represents the 
effect of granulation. The main purpose of this is to prevent elements from 
becoming stationary once they reach the network between supergranules, as our 
supergranular flow proflle does not currently evolve in time. 
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We now update the positions of every element according to its velocity: 

X = X + v^dt, 

y = y + vyrff. 

The next stage is to check which elements have just started or finished emerg- 
ing, cancelling, coalescing or fragmenting. Elements are added and removed 
from the simulation as necessary, and parameters updated. We also check that 
no elements have moved out of the computational domain. If they have, their 
coordinates are translated so that they reappear at the opposite side of the box 
due to the periodic boundaries. 

To construct synthetic magnetograms for the line-of-sight component of the 
magnetic field, we recompute analytically at each time step from the infor- 
mation stored for each magnetic element. This simulated magnetogram is output 
to a file every time step. The method is as described in Section [Ol 

3. Results 

Using the techniques described above, nine simulations are run keeping all in- 
teraction and evolution parameters fixed at values determined to produce the 
most realistic evolution. The only parameter that is varied is the lowest flux 
value taken by emerging bipoles. For the model, we consider a quiet region of 
the solar surface which has no contribution from active regions. The parameters 
common to each model are given in Table [T] The locations of supergranular cells 
and the parameters sud Rq that specify their strengths are also kept fixed, 
as described in Section [2.2l In the nine simulations, the varying values taken for 
the minimum absolute flux of our distribution of emerging bipoles are: 

0min = 10^^ X [4, 8, 10, 20, 30, 40, 60, 80, 100] Mx. 

The results of the magnetic carpet simulations and the effect of varying ^jj^in 
are considered in the following four sections. Two movies showing the full flcld 
of view of the synthetic magnetograms accompany this paper. The movies show 
the = 8 X 10^^ Mx ('mag4.mpg') and ^min = lO" Mx ('magl.mpg') 

simulations between the 50th and 60th hour. By this time, both simulations 
have reached an equilibrium state in which the rate of emergence roughly equals 
the rate of cancellation of magnetic flux. Section 13.11 shows some still images 
from the simulation with (f>^[jj = 8 x 10^^ Mx. 

3.1. Example Synthetic Magnetograms 

Figure IHl shows contour plots of for a simulation in which newly emerging 
bipoles have a total flux in the range 8 x 10^^ - 1 x 10^^ Mx. The images are 
taken at t = 0.1, t = 10, t = 20, t = 225, t = 230 and t = 249 h. These images 
may be compared with the flow profile in Figure IHJb) to compare the location of 
the magnetic elements with respect to the underlying supergranule fiow pattern. 
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Figure 9. Synthetic magnctograms for a simulation with a flux emergence range of 8 X 10^^ 
- 1 X 10^^ Mx. In each frame, red contours represent positive magnetic field, blue contours 
represent negative magnetic field, where 10 contour levels are shown for each polarity with an 
absolute peak value of 70 G. The time in hours at which each image is taken and number of 
individual elements (ns) composing the magnetogram are as follows: (a) t = 0.1, ns=102, (b) 
t = 10, ns=550, (c) t = 20, ns=607, (d) t = 225, ns=595, (e) t = 230, ns=565 and (f) t = 249, 
ns=603. 
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Table 1. Key parameters within the magnetic carpet model, along with their values 
specified. 



Parameter 


Value 


Description 




0, 50 Mm 


x-range for computational box. 


Urmn' J'max 


0, 50 Mm 


y-range for computational box. 




15 000 


Length of simulation. 


stepjength 


1 min 


Number of mins 1 time step represents. 


0max 


10^^ Mx 


Maximum absolute flux for newly 






emerging bipolcs. 


</>o 


lOi*^ Mx 


Minimum unit of flux. 


^rad 


1.5 


Dcflnes distance a magnetic element will travel 




from its initial location under emergence. 


'^rad 


0.5 


Dcflnes interaction distance for cancellation 




and coalescence. 


■^rad 


1.5 


Defines distance a magnetic element will travel 




from its initial location under fragmentation. 


vsg 


0.5 km s~^ 


Peak value for supergranular fiow profile. 




1.5 X 10"'' s-l 


Fragmentation rate. 




8 X 10^** Mx 


Fragmentation parameter. 



At t = the box is empty of magnetic flux. As the simulation progresses, 
bipoles are allowed to emerge following a supergranular cell pattern. After 6 
mins (image (a)) small magnetic elements have begun to appear within the 
supergranular cells. At (b) t = 10 h, many of the elements have begun to en- 
counter one another and interact at the cell edges, and by (c) t = 20 h a magnetic 
network of larger elements has formed along the supergranule boundaries. By 
this stage, the model has reached a steady state in which there is a balance 
between the rates of emerging and cancelling flux. As a result, a similar network 
of magnetic flux exists in each of the later images, (d), (e) and (f). Since the su- 
pergranular flow profile is steady throughout the simulation, the spatial location 
of the magnetic network does not vary greatly; however the exact distribution of 
magnetic flux elements significantly changes. From this it can be seen that the 
steady flow profile docs not lead to the formation of unphysically large magnetic 
elements. However, it would not lead to the random walk of magnetic elements 
across the solar surface as proposed by |Leighton| (|1964p . Since we are presently 
only considering a small, localised area, such a random walk is not important. 
Comparing images (d) and (e) a.t t ~ 225 and t = 230 h, it can be seen that the 
distribution of magnetic elements varies significantly between the two images, 
even though they arc only spaced 5 h apart. This implies that the photosphcric 
recycle time of the simulated magnctograms is short, as desired. A plot of total 
absolute flux versus time is given in Figure llOf b) and discussed in Section 13.21 
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3.2. Total Absolute Flux and Mean Magnetic Field 



For all simulations, the computational box is initially empty of magnetic flux. As 
the simulation progresses, new magnetic bipoles emerge within the supergranular 
cells and magnetic elements interact with one another. It is important to verify 
that the model reaches a steady state in which the absolute flux density oscillates 
about an average value. Plots of the total absolute flux within the box against 
time are shown in Figure [TU] for the four lowest values of (pmin- Clearly, the 
lower the value of the larger the range of flux values for emerging bipoles, 

and thus the greater the rate of flux emergence. This in turn leads to the total 
absolute flux for the simulation levelling off at a higher value. 

It can be seen from Equation ([5]) that N{(j)) will be greater for smaller values 
of (j). Therefore a lower value of 4>Ynin results in a large number of small bipoles 
emerging. For the case where (/iiiiin = 4 x 10^^ Mx (Figure [TOl^a)), the absolute 
flux levels off after approximately 24 h. After this it oscillates between 1.58 x 10^° 
Mx and 1.74 x 10^° Mx. This variation is around 5% from the mean value 
of 1.66 X 10^0 Mx. When ^n^i^ = 8 x 10^*^ Mx, the absolute flux oscillates 
slightly more than this, but the total flux still becomes roughly steady after 
24 h. For higher 4'Ynin- ^^"^ hence a smaller flux emergence range, the variation 
of the absolute flux is much more erratic. For these cases, small numbers of 
large elements emerge, which have a significant effect on the value of the overall 
absolute flux within the box at any one instant in time. When ^j^in = 2 x 10^^ 
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Table 2. Mean values for each simulation for the total absolute flux, absolute flux 
density and number of magnetic elements within the box. 



(xl0i6 Mx) 


Mean total 
absolute flux 
(xlQl" Mx) 


Percentage 
variation of total 
absolute flux 


^fcan absolute 
flux density (G) 


Af can number 
of elements 
per frame 


4 


16.58 


4.9 


6.6 


1497 


8 


10.14 


7.4 


4.1 


570 


10 


8.50 


10.7 


3.4 


430 


20 


7.07 


13.4 


2.8 


211 


30 


9.03 


19.1 


3.6 


153 


40 


5.93 


18.4 


2.4 


120 


60 


6.35 


15.2 


2.5 


88 


80 


5.88 


24.7 


2.4 


74 


100 


5.03 


15.8 


2.0 


63 



Mx fFigurefTUTd)) there are several points where the total flux value spikes due 
to a large bipole emerging, or drops due to the lack of emergence of new bipoles 
and the disappearance of two or more large flux elements through cancellation. 
However, it is clear from Figure [TU] that if c^jnin is chosen low enough, a steady 
state is soon found. 

In several of the simulations with (/)jnin larger than 2 x 10^'' Mx, the total 
absolute flux within the box steadily increases throughout the 250 h period, and 
hence a steady state is not reached. The absolute flux for these simulations may 
level off if the simulation were run for longer, but there would still be a large 
scale flux variation about this value. The mean values for the total absolute 
flux averaged over each simulation, and their percentage variation are shown in 
Table [5] The means are computed for values taken after 30 h of the simulation 
have elapsed, so that the model has had a chance to reach a steady state. The 
highest flux values and lowest variation are found for the lowest value of (j^mm- 

Graphs of the mean magnetic field strength (red line) versus time are shown 
in Figure [TT] for three of the simulations: (a) ^jnin = 4 x 10^^ Mx, (b) ^j^in = 
8 X 10^*^ Mx and (c) 0jnin = 1 x 10^^ Mx. Note that in each graph a different 
scale is used for the vertical axis. Over-plotted on each graph is the number of 
magnetic elements per frame (black line), divided by 100 so that the curve is in 
the same range as the mean field. As one would expect, the mean magnetic field 
strength is higher in simulations where emergence occurs over a larger range. 
The time averaged values for the mean magnetic field in each simulation are 
given in Table [5] along with the mean number of magnetic elements. The mean 
values for the simulations with a (/'min < 3x10^^ are within the observed range of 
3 — 10 G (Gauss) ( [Cranmer and van Ballegooijen, 2010[ ). The mean flux densities 
calculated for (p^nm > 4 x 10^^ Mx are slightly lower than observed values, but 
this is to be expected since realistically, small-scale flux emergence on the Sun 
is not restricted to such a limited range. 
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Time (hrs) 
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Figure 11. (a)-{c) Red line: mean magnetic field (G) and black line: total number of magnetic 
elements (7V/100) versus time for (p^i^ of (a) 4 X 10^*^ Mx, (b) 8 X 10^'^ Mx and (c) 1 X 10^** 
Mx. (d) Red fitted line: time averaged mean magnetic field (G) and black fitted line: average 
number of magnetic elements per frame (N/lOO), for each simulation 



The number of magnetic elements that exist within the box per frame reaches 
a steady state very rapidly in all nine simulations, levelling off after just a few 
hours. This occurs shortly after the first magnetic elements within the simulation 
reach the network between supergranules and begin to coalesce and cancel. The 
average supergranular crossing time in the model, the time taken for a magnetic 
element to reach the edge of the supergranule once it has emerged, is around 3-4 
h. A typical distance that a bipolc might emerge from the edge of a supergranule 
is 7 Mm, if we assume an average velocity of 0.5 km s~^ then each polarity would 
reach the boundary after 3.9 h. The levelling off of the total absolute flux and 
number of magnetic elements is therefore determined by the time scale of the 
flow profile. 

Figure [TlT d') shows the mean field and number of elements (/lOO) averaged 
over the simulation versus (t>xain- Initially, both the mean field and number of 
elements drop very rapidly with increasing (t>Ynin^ then level off around ipjxnn ~ 
2 X 10^^ Mx. 



3.3. Rates and Frequencies 



Within this paper, the rates of magnetic flux emergence and cancellation are 
defined in terms of Mx cm~^ day~^. The frequency of magnetic flux emergence 
and cancellation is defined to be the number of occurrences in cm~^ day"-'^. 
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Figure 12. (a) Emergence (solid line) and cancellation (dashed line) rates for the simulation 
with <^inin = 8 X lO"*^^ Mx. (b) Averaged emergence (black) and cancellation (red) rates for 
all nine simulations. The purple dashed line represents the emergence rates as calculated by 
Equation (O. (c) Emergence (solid line) and cancellation (dashed line) frequencies for the 
simulation with (p^in = 8 X 10^® Mx. (d) Ratio of average emergence to average cancellation 
frequencies for each simulation. 



Figures [T^ a) shows a plot of the variation of the rates of magnetic flux 
emergence (solid line) and magnetic flux cancellation (dashed line) during the 
'f'min = 8 X 10^^ Mx simulation, in Mx cm~^ day~^. In all cases we find that the 
emergence and cancellation rates quickly become approximately equal, confirm- 
ing that the model has reached a steady state. When (/>min = 8 x 10^^ Mx or 
less, the plots are steady, as a large number of small magnetic elements emerge 
and subsequently cancel, helping to keep the emergence and cancellation rates 
steady. As (/)jjiin increases these curves become less steady, but the time averaged 
emergence and cancellation rates for each simulation are all of similar values. 

In simulations with higher ^jnin one would expect the average emergence rate 
to be slightly higher than the average cancellation rate, since the total absolute 
flux in several of these cases increases slightly throughout the simulation. This 
can be seen in Table [31 in which the average emergence and cancellation rates 
for each simulation are listed. Figure [HKb) shows a plot of the average emer- 
gence and cancellation rates for each simulation. The analytical solution for the 
emergence rate calculated from Equation ([9|) is over-plotted in purple. 

Figure [T^] (c) shows a plot of the variation of the emergence and cancellation 
frequencies throughout the simulation for (^min = 8 x 10^^ Mx. The solid line 
represents emergence, the dashed line represents cancellation. It can be seen 



SOLA: km_paperl.tex; 5 August 2011; 0:09; p. 22 



Simulation of Synthetic Magnctograms 



Table 3. Time averaged emergence and cancellation rates and frequencies for each 
simulation. 





Flux 


Flux 


ErQcrgciicc 


C&nccll&tion 




emergence 


cancellation 


frccjUGncy 


frccjUGiicy 


(xlQis Mx) 


rate (Mx cm"-^ 


rate (Mx cm^^ 


(xlO-i« 


(xlO-18 




day-i) 


day-i) 


cm^^ day"-"^) 


cm^^ day^-"^) 


4 


107.37 


107.30 


1503.9 


1540.6 


8 


55.52 


55.41 


344.2 


406.0 


10 


45.69 


45.59 


222.1 


279.7 


20 


25.66 


25.52 


60.4 


98.3 


30 


18.56 


18.32 


29.1 


56.9 


40 


14.61 


14.50 


17.4 


39.6 


60 


10.34 


10.27 


8.5 


23.2 


80 


8.00 


7.89 


5.1 


16.2 


100 


6.53 


6.45 


3.4 


12.5 



that cancellation events occur more often than emergence events even though 
the two processes have similar rates. This is because large flux elements may 
fragment, which leads to more elements that may cancel. As with emergence 
and cancellation rates, the curves are steadier when 4>^\^ is lower. The average 
emergence and cancellation frequencies for each simulation are given in Table |31 
In agreement with Parncll, (.2001|) , the cancellation frequency is always greater 
than the emergence frequency. 

The difference between the frequencies of emergence and cancellation in- 
creases with increasing ipxaxn- This is particularly apparent for larger values 
of ipmim where emergence only produces a relatively small number of large 
elements. These elements may then break apart multiple times to produce large 
numbers of small continually cancelling and coalescing elements, resulting in 
a significantly higher frequency of cancellation than emergence. This occurs 
because emergence within each simulation is restricted to a narrow range of 
flux values, whereas cancellation may occur on any scale. Figure I12f d) shows 
the ratio of emergence to cancellation frequencies. Note that even though the 
cancellation frequency (in cm~^ day~^) may be much greater than the emergence 
frequency, the rates of emergence and cancellation (in Mx cm~^ day~^) are still 
roughly equal. The same quantity of flux emerges and disappears throughout the 
simulation, only the size of the magnetic elements involved in these processes 
changes. 

3.4. Distribution of Flux 

Figures fTSl a') and (b) show number density plots for the number of magnetic 
elements as a function of their flux for the simulations with (j)Ynm = 4 x 10^® Mx 
and ^jnin = 8 x 10^® Mx. In each case the black line represents the number density 
of magnetic elements versus their absolute flux value. This curve is produced by 
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Figure 13. Number densities for (a) flimin = 4 X lO^® Mx and (b) (p^i^ = 8 X 10^'' Mx 
simulations; the black line represents the number density of flux elements versus absolute flux 
content. The blue line illustrates number of emergence events versus absolute flux emerging in 
an event. The red line illustrates number of cancellation events versus absolute flux lost in a 
cancellation event. The green line illustrates number of coalescence events versus absolute flux 
in a coalescence. The purple line illustrates number of fragmentation events versus absolute 
flux of the original element that fragmented. The black dashed line indicates the value of 4>mia- 



analysing the distribution every 600 time steps, which is equivalent to every 10 h, 
and summing up all the flux distributions taken for a particular simulation. The 
length of time between samples was chosen so that the flux within the model will 
have recycled from one sample to the next. This prevents elements from being 
counted more than once. 

The other four lines on each plot represent the occurrence of each of the four 
flux evolution processes described in Section [21 In the simulations it is possible 
to keep track of every event that occurs, so the sample size is a lot larger than 
the flux distribution sample size, as any one element may undergo a number of 
processes. The blue line corresponds to the number density of emergences versus 
the total absolute flux of the emerging bipole, while the red line represents 
the number density of cancellations versus flux removed during the cancellation 
event. It can be seen that the number density of cancellation events is greater 
than the number density of emergence events until (j> > 10"'^^ Mx. This fits with 
previous results that show many more cancellation events occur than emergence 
events, where the events tend to occur between smaller magnetic elements. Above 
a flux of 10^* Mx. fragmentation (purple) is so strong that it affects the elements 
before they can cancel (red). Since the rates of emergence and cancellation are 
roughly equal in each simulation, more large elements must emerge than cancel. 

The green line represents the number density of coalescence events versus the 
absolute flux of the two coalescing elements. In all nine simulations (only two are 
illustrated here) this curve appears to follow the flux distribution curve (black) 
very closely, although the fit is not so good at higher (j) (larger elements do not 
exist for very long within the simulation, so many are missed when sampling the 
flux distribution only once every 600 time steps). This supports the suggestion 
of [Thornton and Parnelll (|201ip . that coalescence is the dominant process for 
small-scale elements. They come to this conclusion because this is one of the 
reasons why fewer small-scale elements are observed on the solar photosphere 
than are found to emerge. It can be seen from the plots in Figures llSf a) and 
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Figure 14. (a) For the simulation with tpmin of 4 X 10^® Mx. The black line represents flux 
density: log plot of frequency of occurrence versus absolute flux of magnetic element. The 
red line represents the fitted line for the full range of flux values, slope=-2.68. The blue line 
represents the fitted line for the range (/> = [10^^, lO'^^] Mx, slope=-1.85. (d) Slopes for flux 
distribution in the range (j> = [10^^, 10^*] Mx for each simulation. 



(b) that the number density of emerging events (bhie) is greater than the flux 
distribution (black) for small values of 0, supporting the fact that coalescence 
is dominant for small-scale fields. In order for our model to be more conclusive 
on this matter however, it is necessary to allow the existence and emergence of 
even smaller flux elements than 0o = 10^^ Mx. Since the green (coalescence) 
line follows the flux distribution (black) so closely, this implies that any and all 
magnetic elements may undergo coalescence, it is not flux dependent. 

Fragmentation is defined in our simulations to be strongly flux dependent. 
The purple line represents the number density of fragmentation events versus 
the flux of the original magnetic element. The curve for fragmentation is lower 
than those of all of the other processes for small (j>. However, it can be seen 
that as an element's absolute flux increases the fragmentation process becomes 
dominant. 

Figure iMT a) is a plot of the frequency of occurrence of flux elements versus 
their flux content, for the simulation with ^min = 4 x 10^^ (black line). Units 
of xIO^"'^ Mx^^ cm~^ have been chosen for the y-axis in order to compare the 
plot with Figure 5 of IParnell et al. \ (|2009p . Two straight lines have been fitted 
to the data. We note that our data do not span three orders of magnitude, as 
is technically required to compute a power law relation. However, power laws 
have been used for observational data so we have fitted one for comparison. The 
red line is fitted using all of the data, and has a slope of -2.68. The blue line 
is fitted for elements in the range 10^^ — 10^^ Mx, therefore missing the low 
emergence upper flux range, as the sample size for our model is much better for 
smaller elements. Note that as we are considering an isolated region of the solar 
surface, the upper bound is solely dependent on emergence and coalescence. No 
elements may enter the domain due to the dispersal of magnetic flux from active 
regions. As a result, we find that our distribution of elements falls off faster than 
is observed at higher values of fiux. Hence we restrict the power law index to 
those flux elements that we adequately model. The fit of the blue line is good 
in this range. It has a slope of -1.85, which is in agreement with the findings of 
IParnell et aLI(l2009p . They fit a line with a slope of -1.85 to data for fiux elements 
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observed by SOHO/MDI and Hinode /SOT that spans the range 10^^ - lO^^ Mx, 
therefore our model, which only models the lower range of this, fits this well. The 
method that they use to detect magnetic fiux features also differs from the way 
that we define them. They use a 'clumping' method, whereas our definition of 
a magnetic element is more suited to the 'downhill' method of feature tracking 
(De Forest et a/. IMTTl Lamb et al. . nUMHUWl . The downhill method is good at 
picking out individual peaks in magnetic flux. Particularly for larger magnetic 
regions, several peaks may be detected by the downhill method and counted as 
separate magnetic features. However, the same region may be counted as a single 
feature by the clumping method (see Figures 1 and 7 of lDeForest et aLI(|2007p ). 
A flux distribution detected using downhill would result in a 'tail' towards higher 
fluxes, as is seen in our Figure [Hla). A feature tracking study of our synthetic 
magnetograms would be of interest, using a clumping method such as is described 
bv lThornton and Parnelll (POTT|) . 

The equation of IParnell et al. \ (|2009p describing the frequency of occurrence 
of elements with an absolute flux of has the form 



where Nf = 3.6 x 10~^ is the value obtained using a clumping method of feature 
tracking^. In our simulations, we obtain iVf = 1.2 x 10~^^, which is 3.3 times 
larger. A possible reason for this is the difference between our definition of a flux 
element, and the way that features are deflned bv IParnell et aLI(|2009p . As can 
be seen in their Figure 4(a), the clumping method they use produces relatively 
large, irregularly shaped 'flux massifs'. Since we define magnetic elements to be 
compact and circular, such a feature in our model would be composed of many 
elements. For this reason we require many more small magnetic elements to 
describe the flux distribution of the quiet Sun, hence the larger value of Nf in the 
number density equation. It should be noted that while our definition of a mag- 
netic element is different from that of IParnell et al.\ (|2009p after it has emerged, 
this does not affect the emergence power law ( [Thornton and Parnell, 20 lip , as 
during the process of emergence our definitions of magnetic elements agree. 

Our reasoning for the difference in values of is supported by some simple 
calculations. If a = 1.85, then the average fiux density (i^avg) and average fiux 
of a magnetic element ((/)avg) are given by 




(13) 




(14) 



and 



'avg - 




(15) 



■^C. E. Parnell (private communication). 
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Taking our values of TVf = 1.2 x 10"^'' cm~^, (p^^^ = 4 x 10^^ Mx and (/)max = 
10^* Mx we find 0avg = 1-5 x 10^'^ Mx and B^vg = 6.11 Mx cm^^ Our B^vg 
is consistent with the values given in Table [51 and is realistic for the quiet 
Sun. If we now consider the iVf = 3.6 x 10"" cm'^ of IParnell eTflD ((20091) . 
but limit the range of flux values considered to be consistent with the quiet 
Sun (4 X IQis - 1020 Mx), we get (j)avg = 5.0 x lO^^ Mx and B^vg = 6.6 Mx 
cm~2. Our absolute flux density is very similar to theirs, but our average flux 
element is much smaller. We also find that we have around three times more 
magnetic elements per unit area than IParnell et aT1(|2009p (iVtot = 4.06 x 10~^^ 
cm~2 compared with iVtot = 1-30 x 10~^^ cm"^). Therefore the difference in 
the parameters of the power law is due to our varying definition of 'magnetic 
elements' to 'magnetic features'. 

The slope of the fitted line in the range 10^^ — 10^^ Mx for each simulation 
is plotted versus 0jnin in Figure IHTb) . One can see that as (f>^^i^ increases the 
slope becomes less steep. This is because as only larger elements emerge, these 
large elements have more of an impact on the number density. Such a limited 
range for emerging bipoles is less realistic, so we would expect this to be less of a 
match with the results obtained through observations. When the lower bound of 
emerging fiux is less than 10^^ Mx, the power law index converges around -1.8. 

3.5. Lifetime of Magnetic Elements 

The plots in Figure [15] relate to the lifespan of magnetic elements within each 
simulation. Within our model, an element is defined to 'die' when its flux changes, 
either by fragmenting into two new elements or cancelling or coalescing with 
another element. An element begins its life when it newly emerges; has just 
split from another as a result of fragmentation; or is produced by two separate 
elements cancelling or coalescing together. 

Our definition of an element's lifespan likely differs somewhat from that of 
an observer studying magnetogram data, as does our counting of magnetic ele- 
ments. It is easy to count and keep track of the processes that magnetic elements 
undergo within our model because they are treated as individual point sources. 
However, when the synthetic magnetogram is created using the method described 
in Section I^TTl many of these elements overlap to produce fewer, larger magnetic 
elements. In addition, the resolution and cadence of real data are not always high 
enough to be able to detect the smallest and fastest evolving elements. This is 
another reason why further study of the photospheric model using the synthetic 
magnetogram series as input into the variety of feature tracking techniques that 
have been produced would be of interest. 

Figure [TST a) shows the log frequency of occurrence of magnetic elements with 
a lifespan of t h, for the simulation in which (^juin = 4 x 10^^ Mx. The mean 
lifespan of all elements is plotted for each simulation in Figure fTST b). It can 
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Figure 15. (a) Frequency of occurrence of magnetic elements with a lifespan of t hours, 
within the (/<niin = 4 X 10^^ Mx simulation, (b) Mean lifespan of a magnetic element in hours 
for each of the simulations, (c) The age of the oldest existing magnetic element throughout 
the simulation, for flimin = 4 X 10^^ Mx. The red dashed line represents the mean age of the 
oldest element averaged over the simulation, (d) Black stars: mean maximum age of an element 
averaged over each simulation. Blue line: mean maximum age of an element averaged over all 
simulations. Red stars: mean lifespan of a magnetic element within each simulation. 

be seen that most elements do not even 'live' for one hour before fragment- 
ing, cancelling or coalescing with another. The mean lifespan of an element is 
only around 9-20 mins. Observationally, it makes more sense that larger, more 
isolated magnetic elements would be long-lived. However, since fragmentation 
is so frequent within our model, this is not the case for our simulations. The 
longest lived elements most likely occur at the beginning of the simulation before 
the number of magnetic elements increases to the point where they interact 
frequently; or are isolated small elements that do not fragment because their 
flux is so low. Within the synthetic magnetograms, large, irregularly shaped 
features form where several individual elements lie close together but have not yet 
coalesced. These elements tend to appear in the network between supergranules. 
Examples of these can be seen in the six images shown in Figure [SI 

Figure fTSTc) shows a plot of the age of the oldest magnetic element existing 
within the model versus time, within the 4>^[^ = 4 x 10^^ Mx simulation. The 
dashed red line represents the mean maximum age averaged over the simulation. 
The large spikes represent occasional long-lived magnetic elements, but in general 
the maximum age of an element throughout the simulation remains at roughly 
3-4 h. In every case the mean maximum age computed for the whole simulation 
is around this value, as can be seen in FigurejTSJd). This is similar to the lifetime 
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Figure 16. A simulation with emergence of bipoles in the range 8 X 10^® — 10^^ Mx. Emergence 
is switched off after 50 h (3 000 time steps), and the simulation is allowed to run for a further 
50 h. (a) Total absolute flux within the simulation versus time, (b) Mean field and number of 
magnetic elements (A'^/100) versus time, in red and black respectively, (c) and (d) show the 
synthetic magnctograms produced at t = 50 h and t = 60 h. 



of ephemeral regions, determined by [Harvey and Martin (jl973p to be around 4.4 
h. I Zhou et al. \ (j2010[) find the hfespan of internetwork features to be between 1 
and 20 mins. with a mean of 2.9±2.0 mins. The mean Hfespan for magnetic 
elements within our model is 9-20 mins. 

The photospheric recycle time is calculated by dividing the mean field by the 
emergence rate. For our two most realistic simulations, 0niin = 4 x 10^^ Mx and 
'?^min = 8 X 10^^ Mx, we find the recycle time to be 1.48 h and 1.75 h respectively. 
This is in excellent agreement with |Hagenaar et al.\ (|2008p 's recycle time of 1-2 
h. 

The next section considers one of many possibilities for future studies using 
this theoretical model for the magnetic carpet, in this example, emergence is 
switched off midway through the simulation. 
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3.6. Switching Off Emergence 

Clearly, many other studies may be undertaken using our theoretical model 
simply by varying different input parameters to determine their effect upon 
photospheric evolution. One interesting test is to 'switch off' emergence at some 
stage in the simulation to observe how rapidly flux disappears. An example of 
such a simulation was run with identical parameters to those described in Sec- 
tion[3l The flux emergence range for newly appearing bipoles is 8 x f 0^^ — f x f 0^^ 
Mx. Emergence is switched off at t = 50 h, and the simulation is allowed to run 
for a further 50 h. Results from this experiment are shown in Figure [TOl Plot (a) 
shows the total absolute flux and (b) shows both the mean field and the number 
of elements (N/lOO) throughout the simulation. It can be seen that all three of 
these quantities rapidly decrease as soon as emergence is switched off after 50 
h, levelling off and becoming steady at roughly t = 60 h. The total number of 
magnetic elements becomes very small as all remaining flux within the computa- 
tional box rapidly cancels and coalesces together at the supcrgranulc boundaries. 
The synthetic magnetogram images (c) and (d) arc taken at t = 50 and t ~ 60 
h. In image (d) one can see the few remaining magnetic elements. They are 
spaced far from one another where the supergranular velocities are small at the 
boundaries, so they do not encounter one another before the simulation ends. 
If the supergranular flow profile was also allowed to evolve in time, this would 
likely aid the process by increasing and randomising further the motion of the 
remaining magnetic elements, eventually removing them completely from the 
simulation. 



4. Discussion, Conclusions and Future Work 

The aim of this paper is to construct a realistic model for the photospheric 
evolution of the solar magnetic carpet, which we intend to use as a lower bound- 
ary condition of 3D non-potential modelling. We have built into this model the 
processes of flux emergence, cancellation, coalescence and fragmentation, as well 
as a steady supergranular fiow profile that infiuences the motion of magnetic 
elements. Many parameters for the model are taken from studies of observa- 
tional solar data, such as the probability distribution for newly emerging bipoles 



Simon and Leighton, 1964 



Thornton and Parnell, 20111 and the peak radial velocity of a supergranule 



'aniveni et al, 2004| . 



A series of nine simulations of length 250 h were run, keeping all parameters 
fixed apart from the minimum value of total absolute fiux, (jj^a^^, that a newly 
emerging bipole may take. The maximum value for newly emerging bipoles was 
fixed at (/>max = 10^^ Mx. The lower the value of (f>^[ji is taken to be, the larger 
the range of emerging flux elements is, and hence the more flux is emerged into 
the system. The upper bound produces a cap on the size of large elements, and 
we assume that we arc considering a quiet area of the Sun, with no input from 
active regions. A larger flux emergence range is more realistic so as one would 
expect, the lower the value of ^min, the closer the simulation results agree with 
solar observations. 
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Particularly for 0niin = 4 x 10^^ Mx and (j)^^^ ^ 8 x 10^^ Mx, the total 
absolute flux within the simulation quickly reaches a steady state in which 
the rates of emergence and cancellation are roughly equal. These simulations 
also result in a mean magnetic field that is within the range determined from 
observations. For the less realistic simulations with the highest (/)niin (^^ 3 x 10^^ 
Mx) the total absolute flux takes longer to reach a steady state illustrated by the 
fact that the mean emergence rate is slightly higher than the mean cancellation 
rate. For low ^mim mean field is 4-6 G, which fits observational data. 

The number of magnetic elements within each simulation very quickly reaches 
a steady state in all cases however. This occurs shortly after the first magnetic 
elements reach the boundaries of the supergranules and begin to interact with 
one another. 

Although the cancellation and emergence rates, defined in Mx cm~^ day"-'^, 
become roughly equal for all simulations, it is also the case that the cancellation 
frequency, in cm~^ day~^, is always greater than the emergence frequency. This is 
in rough agreement with the theoretical magnetic carpet model of lParneIII(|200ip . 
who attribute the difference to large numbers of cancellations between small 
magnetic elements arising through fragmentation. They suggest that the energy 
release from so many small cancellation events could significantly contribute to 
the energy required to heat the solar corona. In future in our 3D simulations, 
we will consider the energy build up at such locations. 

The flux distribution resulting from the simulations that include a larger 
range of fiux emergence are in agreement with the power-law distribution of 
IParnell et aLI(|2009p . The slope of their fitted line is -1.85, while the slopes pro- 
duced by our simulations with ^min = 4 x 10^^ Mx and (f>^[j^ = 8 x 10^^ Mx are 
-1.85 and -1.82 respectively. Our definition of magnetic elements means that our 
fiux distribution would fit better with a downhill method of feature tracking. A 
feature tracking study of our synthetic magnetograms using a clumping method 
would be useful to check the power law, in addition to other values such as 
lifetimes and physical extent of magnetic elements. 

Our model produces a highly dynamic small-scale photosphere as desired, 
with the mean lifespan of a magnetic element in any of the simulations being 
just 9-20 mins. We also find a photospheric recycle time of just 1.48 h, this is in 
good agreement with |Hagenaar et a^.| (|20(]5)) 's recycle time time of 1-2 h. 

As suggested in Section 13.61 an evolving supergranular flow profile is one im- 
provement that could be made to the model. Simon and Leighton (jl964p deter- 



mine the average lifespan of a supergranule to be 20 h, whereas Wang and Zirin 



(|1989p find their lifetime to be > 50 h. ISimon et aL\ (I200ip assign hfetimes of 
between 18 and 42 h to the supergranules within their model. In any case, these 
studies show that a magnetic carpet model that runs for longer than around 
a day should consider a time evolving flow pattern. The flow profile evolution 
could either be built into the model or taken from observational data using a 
method such as described by |Potts, Barrett and Diver] (j2004p . They track solar 
photospheric flows by using a ball tracking technique. 

The fragmentation process is limitation of our model. Rather than allowing 
fragmentation to arise naturally as a consequence of photospheric flows, we 
currently impose it. This is another feature of the model that could be improved 
in future. 
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The next stage in our process of constructing a non-potential coronal model 
for the solar magnetic carpet is the theoretical study of small-scale basic interac- 
tions between magnetic flux elements, using the technique of |van Ballegooijen, Priest and Mackay[ 
(|2000p . This will produce a continuously evolving non- linear force-free magnetic 
field in response to photospheric motions ( [Meyer et al., 2011b I . This method has 



recently been applied to the evolving structure of magnetic filaments ( Mackay and van Ballegooijen, 2009| ,| 



and a decaying active region observed by SOHO/MDI ( |Mackay et al, 2011D 



It is useful to understand these simple interactions before moving onto more 
complex cases. We will later insert the synthetic magnetograms produced by 
this 2D model into the 3D coronal model as lower boundary data. 

We will study many aspects of the resultant coronal magnetic field, such 
as locations of electric currents, free magnetic energy, and coronal null points, 
as these are associated with the problem of coronal heating. It would also be 
of interest to conduct this study side by side with an identical study of SDO 
magnetogram data, for which the coronal magnetic field could also be modelled. 
Since we know exactly which processes are occurring where within the simulated 
synthetic magnetograms illustrated here, we would be able to relate these with 
events occurring in the real magnetogram data, and study the simulated coronal 
field in both cases. 

We conclude that we have successfully produced a realistic model for the 
photospheric evolution of the solar magnetic carpet that reproduces many ob- 
servational properties. This will now be used to produce a 3D model for the 
small-scale coronal magnetic field under controlled circumstances. 
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